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HIGHLIGHTS 


•  SSM50  has  better  oxygen  storage  capacity  than  LSM50. 

•  The  absorbed  02  species  on  SSM50  (100)  surface  are  readily  dissociated. 

•  SSM50  has  the  best  catalysis  for  the  ORR  among  the  considered  systems. 
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Cubic  perovskite  SmxSri_xMn03  (SSM)  surface  and  bulk  models  have  been  constructed  to  simulate  the 
oxygen  reduction  reactions  by  employing  first-principles  calculations.  The  results  demonstrate  that 
oxygen  vacancies  can  be  formed  easily  in  Sm0.5Sr0.5MnO3  (SSM50).  The  oxygen  migration  barrier  in  bulk 
SSM50,  which  is  predicted  by  the  nudged  elastic  band  (NEB)  method,  is  the  lowest,  while  the  adsorption 
energy  of  02  molecular  on  SSM50  (100)  surface  is  the  lowest  among  the  considered  doping  systems, 
indicating  the  potential  application  of  SSM50  as  a  cathode  for  intermediate-temperature  solid  oxide  fuel 
cell  (IT-SOFC).  The  reaction  mechanisms  of  oxygen  reduction  on  SSM50  (100)  surface  have  also  been 
studied. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cell  (SOFC)  has  received  increasing  attention  in 
recent  decades  because  it  has  higher  efficiency  and  lower  pollutant 
emission  compared  with  other  conventional  energy  conversion 
devices  [1—3].  Lanthanum  strontium  manganite  (LSM)  is  the  most 
conventional  cathode  for  high  temperature  SOFC  [4-6].  However, 
with  the  decrease  of  operating  temperature,  the  polarization 
resistance  increases  significantly  [7],  which  limits  the  application  of 
the  LSM  based  cathodes  for  intermediate-temperature  SOFC  (IT- 
SOFC).  The  polarization  resistance  of  the  cathode  originates  from 
the  oxygen  reduction  reaction  (ORR)  [8].  The  oxygen  adsorption  on 
the  cathode  surface  and  the  oxygen  vacancy  migration  in  the 
cathode  are  the  crucial  factors  which  affect  the  polarization  resis¬ 
tance  during  the  ORR  9].  These  two  factors  can  be  characterized  by 
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oxygen  adsorption  energy  and  ionic  diffusion  activation  energy, 
respectively. 

Great  efforts  [10-13]  have  been  made  to  study  the  mechanism 
of  ORR  and  develop  SOFC  cathodes  with  good  performance.  Among 
these  studies,  Density  Functional  Theory  (DFT)  calculation  has 
become  an  important  and  powerful  tool  which  allows  us  to  gain 
further  insight  into  the  nature  of  the  oxygen  exchange  reaction. 
Oxygen  adsorption  and  diffusion  in  the  common  used  LSM  cathode 
have  been  studied  systemically  by  Kuklja  [13].  It  is  found  that 
the  ORR  mainly  occurs  on  the  Mn02-terminated  surface  and 
the  surface  migration  of  the  adsorbed  O  species  is  slower  than 
that  of  surface  O  vacancies.  Liu  et  al.  [8,14,15]  have  examined  the 
electrochemical  reaction  activities  of  cathode  with  different 
perovs kite-type  materials  such  as  PrxSri^xMn03,  NdxSri_xMn03, 
PmxSri_xMn03,  and  SmxSri_xMn03  (SSM).  It  is  found  that  only 
(SSM)  exhibits  superior  electrochemical  performance  [14].  SSM  is 
an  electronic  conductor,  but  it  can  become  a  mixed  ionic-electronic 
conductor  (MIEC)  under  cathodic  polarization  because  polarization 
promotes  an  in-situ  electrochemical  reduction  of  the  Mn4+  in  SSM 
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to  Mn3+  with  the  introduction  of  oxygen  vacancies  into  SSM.  The 
previous  work  only  considered  the  case  of  x  =  0.5,  and  there  has 
been  no  discussion  on  the  mechanism  of  oxygen  reduction. 
Therefore,  the  DFT  calculations  on  the  complex  SSM  are  conducted 
in  the  present  work.  Besides,  the  formation  of  bulk  oxygen  va¬ 
cancies  and  their  diffusions,  as  well  as  the  oxygen  adsorption 
property  on  the  SSM  surface  are  studied. 

2.  Computational 

The  ab  initio  DFT  computer  code  VASP  [16,17],  which  uses  the 
plane-wave  pseudopotential  total  energy  calculation  method  [18],  is 
introduced  to  conduct  all  the  first-principles  calculations.  The  ex¬ 
change  correlation  of  electrons  is  parameterized  by  Perdew-Burke— 
Ernzerhof  (PBE)  [19]  formula  within  the  generalized  gradient 
approximation  (GGA).  The  Brillouin  zone  sampling  4x4x4  Mon- 
khorst-Pack  k-points  [20]  for  bulk  and  2  x  3  x  1  for  surface  is 
adopted,  and  the  cutoff  energy  is  set  to  be  400  eV,  according  to  the 
result  of  convergence  test.  Using  the  block  Davidson  scheme,  both  the 
atomic  positions  and  cell  parameters  are  optimized  until  the  residual 
force  is  below  0.01  eV  A-1.  Since  the  LaMn03-based  cathode  material 
has  a  cubic  structure  under  SOFC  operating  conditions  [15,21],  and 
SrMn03  also  has  a  cubic  structure  at  high  temperatures,  we  first 
constructed  an  ideal  cubic  perovskite  structure  of  SrMn03  (space 
group  Pm3m).  The  calculated  lattice  parameter  after  structure  opti¬ 
mization  is  3.802  A,  which  is  in  good  agreement  with  the  experi¬ 
mental  values  of  3.805  A  [22  .  When  part  of  or  all  the  Sr  atoms  have 
been  replaced  by  Sm  atoms,  there  are  only  tiny  changes  of  the  lattice 
parameter  but  the  structures  remain  cubic.  Based  on  the  above  re¬ 
sults,  it  is  assumed  that  SSM  cathode  also  has  a  cubic  perovskite 
structure  in  our  calculations.  Therefore,  only  the  highly  symmetric 
cubic  perovskite  structure  is  considered  similar  to  previous  study 
[14].  All  the  doped  systems  are  constructed  from  a  relaxed  (2x2x1) 
SrMn03  supercell,  with  Sm  doping  in  SSM  atx  =  0,  0.25,  0.5,  0.75, 1, 
respectively.  In  the  calculation,  the  4f  electrons  of  Sm  are  treated  as 
core  electrons,  and  the  electronic  configurations  are  considered, 
included  Sm  (5s25p65d16s2),  Sr  (4s24p65s2),  Mn  (3d54s2),  and  O 
(2s22p4). 

3.  Results  and  discussion 

3.1.  Bulk  diffusion  ofO  vacancy 

A  good  SOFC  cathode  should  be  a  good  MIEC.  The  bulk  diffusion 
of  O  vacancy  in  SSM  is  studied  in  the  first.  The  diffusion  activation 
energy  Ea  is  applied  to  characterize  the  difficulty  of  an  oxygen  ion 
(or  vacancy)  migration  in  bulk  SSM.  According  to  Liu  et  al.  [13],  the 
Ea  is  divided  into  two  parts:  the  oxygen-vacancy  formation  energy 
(Evac),  and  the  migration  barrier  (Em)  of  the  oxygen  ion  (or  vacancy). 
And  they  are  described  as 

Ea  +  Eyac  +  Em .  ( 1 ) 

An  O  vacancy  is  introduced  into  the  SSM  (2x2x1)  supercell, 
forming  a  defect  structure  of  SmxSri_xMn02.75.  The  corresponding 
oxygen-vacancy  formation  energy  can  be  calculated  by 

Evac  =  EsmxSr^MnCb.ys  +2^°2  -  ESmxSr^  xMn03 >  (2) 

where  ESmxsri_xMno2  75  and  ESmxsrl  xMno3  represent  the  total  energies 
of  bulk  SSM  with  and  without  an  oxygen  vacancy,  respectively. 
The  total  energy  of  O2  ( E0l )  is  calculated  in  a  9  A  x  10  A  x  11  A 
box,  to  prevent  the  effect  of  the  degeneracy.  After  geometry 
optimization,  all  the  reduced  structures  show  distortions  around 
the  vacancy  site.  The  calculated  Evac  are  listed  in  Table  1.  It  is  found 


Table  1 

Calculated  energies  for  oxygen-vacancy  formation  Evac  and  the  migration  barriers 
(£m)  in  varies  bulk  SSM. 


SSM 

Ev  ac  (ev) 

(eV) 

Ea  (eV) 

SrMn03 

1.52 

0.73 

2.25 

Sm0.25Sr0.75MnO3 

2.30 

0.85 

3.15 

Sm0.5Sr0.5MnO3 

2.58 

0.75 

3.33 

Smo.75Sro.25Mn03 

3.40 

0.96 

4.36 

SmMn03 

3.75 

1.02 

4.77 

that  the  oxygen-vacancy  formation  energy  depends  strongly  on 
the  dopant  concentration.  As  we  can  see,  Evac  becomes  smaller  as 
the  Sm  doping  concentration  increases.  When  half  of  the  Sr  atoms 
are  replaced  by  Sm,  the  value  of  Evac  is  2.58  eV,  which  is  smaller 
than  that  of  Lao.5Sr0.5Mn03  (LSM50)  (3.25  eV)  [14],  indicating  the 
oxygen  vacancies  can  be  formed  more  easily  in  SSM50  than  in 
LSM50.  In  other  words,  SSM50  has  better  oxygen  storage  capacity 
than  LSM50. 

To  predict  the  migration  barrier  of  an  oxygen  ion  in 
SmxSri_xMn02.75,  we  assume  that  the  oxygen  ion  is  diffused 
directly  to  its  neighbor  vacancy  site,  as  shown  in  Fig.  1.  Under  this 
assumption,  a  migrating  oxygen  ion  has  to  pass  through  the 
triangular  orifice  consisting  of  two  large  Sr  (or  Sm)  and  one  small 
Mn  atoms  in  the  crystal  lattice.  We  first  find  out  the  transition 
states  using  the  nudged  elastic  band  (NEB)  method  [23,24  .  The 
energy  difference  between  Transition  State  (TS)  and  the  Initial 
State  (IS)  represents  the  migration  barrier.  The  migration  barriers 
in  various  SSM  systems  are  shown  in  Table  1,  together  with  the 
corresponding  activation  energy.  It  can  be  concluded  that  the 
oxygen  ion  overcomes  a  barrier  of  0.75  eV  when  migrating  to  its 
neighbor  site  in  SSM50.  The  same  as  the  case  of  other  cathode 
materials,  e.g.  (Ba,  Sr)(Co,  Fe)03  [25],  the  value  of  the  migration 
barriers  varies  with  the  different  doping  concentration. 

According  to  the  Equation  (1 ),  the  diffusion  activation  energy  Ea  is 
the  sum  of  Ev ac  and  Em.  According  to  Table  1,  the  Ea  of  SSM50  is 
3.33  eV,  a  little  bit  higher  than  that  of  SrMn03  and  Sm0.25Sr0.75MnO3. 
However,  because  the  valence  of  Sr  is  +2  in  SSM50,  which  is  lower 
than  that  of  Sm  (+3),  half  of  the  Mn3+  ions  will  change  to  Mn4+  for 
the  balance  of  the  charge.  The  coexistence  of  Mn3+  and  Mn4+  will 
form  an  electronic  transmission  channel  (Mn3+-0-Mn4+),  which 
can  lead  to  a  better  electronic  conductivity  [26]  for  SSM50  than 
SrMn03  and  Sm0.25Sr0.75MnO3  (SSM50  has  more  conducting  elec¬ 
trons  than  Smo.25Sro.75Mn03  does).  Based  on  the  above  consider¬ 
ations,  SSM50  is  most  suitable  as  SOFC  cathode  among  the 
considering  systems. 


£  Sr(Sm) 

J  Mn 

•  O 

O  O  vacancy 


Fig.  1.  Schematic  representations  of  oxygen  ion  migration  in  defect  bulk  SSM. 
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3.2.  Surface  adsorption  and  diffusion 

The  oxygen  adsorption  at  the  cathode  surface  is  a  critical  step 
during  the  cathode  oxygen  reduction  process.  In  order  to  study  the 
adsorption  characteristics,  we  model  the  O2  (and  0)  adsorption  on 
the  SSM  (100)  surface,  since  it  is  energetically  the  most  stable 
surface  among  the  low-index  (100),  (110),  and  (111)  surfaces.  The 
SSM  (100)  surface  consists  of  alternating  MnC^  and  Sr(Sm)0 
plane.  According  to  the  results  of  convergence  test,  a  six-layer  slab 
is  thick  enough  to  simulate  the  adsorption.  Supercells  with  peri¬ 
odicity  (2x1)  are  employed  to  reduce  the  influence  of  the 
interaction  between  the  adsorbed  species,  and  the  periodically 
repeated  slabs  are  separated  by  a  15  A  thick  vacuum  layers  to 
avoid  the  interaction  between  the  two  surfaces.  Fig.  2  shows 
the  side  elevation  of  O2  adsorption  on  SSM  (100)  surface.  All  of  the 
six  layers  are  allowed  to  be  relaxed  during  the  geometry 
optimization. 

The  adsorption  energy  (Eads)»  which  corresponds  to  the  energy 
gained  from  the  adsorption  process,  is  defined  by 

^ads  =  ^adsorbate+slab  —  ^slab  —  ^adsorbate  5 

where  Eadsorbate+siab  is  the  total  energy  of  the  adsorbed  assembly, 
Esiab  and  Eadsorbate  are  the  energies  of  the  clean  surface  and  the 
adsorbate,  respectively.  The  value  of  Eads  is  a  criterion  to  determine 
the  difficulty  and  stability  of  the  adsorption.  The  more  negative  the 
Eads  is,  the  easier  and  the  more  stable  the  adsorbate  can  be  adsor¬ 
bed  on  the  surface.  Previous  study  [13,27]  shows  that  B  cations  are 
more  active  than  A  cations  for  oxygen  reduction  on  several 
perovskite-type  ABO3  cathodes.  Only  the  outermost  Mn  atoms 
(atop  Mn)  on  the  Mn02-terminated  surface  are  considered  in  this 
study. 

The  calculated  adsorption  energy  and  structural  parameters  for 
O2  adsorption  on  the  atop  Mn  of  SSM  (100)  surfaces  are  summa¬ 
rized  in  Table  2.  The  calculated  bond  length  of  O2  molecule  in 
equilibrium  gas  phase  is  1.234  A,  slightly  larger  than  experimental 


J  Sr(Sm) 
J  Mn 

•  O 


Fig.  2.  The  side  elevation  of  02  adsorption  on  SSM  (100)  surface. 


Table  2 

Calculated  adsorption  energy  (Ea ds)  and  structural  parameters3  for  02(0)  adsorption 
on  the  atop  surface  Mn  atoms  of  SSM  (100)  surfaces. 


o2 

0 

Fads 

doi-Mn 

do 1-02 

Fdiss 

Fads 

doi-Mn 

(eV) 

(A) 

(A) 

(eV) 

(eV) 

(A) 

SrMn03 

-3.09 

1.773 

1.275 

0.62 

-2.49 

1.589 

Smo.25Sro.75Mn03 

-3.13 

1.741 

1.281 

0.70 

-2.78 

1.591 

Sm0.5Sro.5Mn03 

-3.94 

1.737 

1.285 

0.53 

-3.35 

1.591 

Sm0.75Sro.25Mn03 

-3.93 

1.728 

1.281 

0.59 

-3.39 

1.592 

SmMn03 

-3.52 

1.730 

1.283 

0.74 

-3.25 

1.587 

doi-Mn:  distance  between  01  and  the  nearest  Mn  atom,  doi-02:  distance  between 
01  and  02. 

a  01  atom  nearest  to  the  surface,  02  atom  far  from  the  surface  (as  shown  in  Fig.  1 ). 


value  of  1.210  A.  When  adsorbed  on  the  SSM  surface,  the  bond 
lengths  of  the  O2  molecule  increase  to  1.275-1.285  A.  The  increase 
of  the  bond  length  designates  that  the  absorbed  O2  species  are 
readily  to  be  dissociated.  The  distance  between  01  and  the  surface 
Mn  ranges  from  1.728  A  to  1.773  A,  which  are  smaller  than  the  sum 
of  the  radius  of  Mn3+/Mn4+  and  O2-,  demonstrating  that  Mn  and  O 
are  connected  by  chemical  bond. 

Besides  the  adsorption  energy,  the  O2  dissociation  energy  is  also 
one  of  key  parameters  that  need  to  be  discussed.  Using  the  same 
method  described  by  Mastrikov  [28],  the  dissociation  energies 
(Ediss)  of  the  adsorbed  02  have  been  calculated  and  listed  in  Table  2. 
The  Ediss  ranges  from  to  0.53-0.74  eV,  indicating  that  the  adsorbed 
O2  can  be  easily  dissociated. 

In  SSM50  system,  the  Eads  and  EdiSS  of  O2  molecular  are  calcu¬ 
lated  to  be  -3.94  V  and  0.53  eV,  lower  compared  to  that  of  any 
other  systems;  and  the  bond  length  of  the  adsorbed  O2  is  larger 
than  that  of  the  other  systems  too,  implying  that  SSM50  may  have 
the  best  catalysis  for  the  ORR  among  these  five  systems.  Along 
with  the  O2  molecule  adsorption,  the  adsorption  properties  for  O 
atom  on  the  SSM  (100)  surfaces  are  calculated.  The  result  shows  a 
strong  preference  for  O  atom  adsorption  on  the  SSM50  (100) 
surface. 

Previous  study  [28]  demonstrates  that  surface  vacancy  diffusion 
energy  (Edifr)  strongly  differs  from  that  in  the  bulk.  Thus,  the  Edifr 
also  has  been  calculated  (Table  3)  using  the  NEB  method.  The  same 
as  the  diffusion  in  the  bulk  SSM,  we  assume  that  the  oxygen  ion  is 
diffused  directly  to  its  neighbor  vacancy  site  (Fig.  3).  It  is  found  that 
Edifr  is  higher  than  the  relevant  Em,  clarifying  that  it  is  surface 
diffusion  that  is  limiting  the  overall  oxygen  vacancy  diffusion  rate. 

Combining  the  oxygen-vacancy  formation  energy,  the  migra¬ 
tion  barrier,  the  adsorption  and  dissociation  energy,  as 
shown  in  Tables  1  and  2,  we  can  conclude  that  SSM50  is  a  po¬ 
tential  SOFC  cathode  material  with  excellent  performance  for 
ORR. 

3.3.  PDOS  analysis  of  the  adsorbed  assemblies 

To  explore  the  deep  mechanism  of  the  interaction  between  the 
adsorbates  and  surfaces,  partial  densities  of  states  (PDOS)  calcula¬ 
tions  on  the  O2  molecular  adsorbed  SSM50  (100)  surface  are  carried 
out.  The  2 p  orbitals  of  the  adsorbed  O  atoms  and  the  surface  O 


Table  3 

Calculated  oxygen  vacancy  diffusion  energy  (Edifr)  on  SSM  (100)  surface. 


System 

SrMn03 

Sm0.25Sro.75 

Mn03 

Smo.5Sr0.5 

Mn03 

Smo.75Sro.25 

Mn03 

SmMn03 

Fdiff(eV) 

0.86 

0.95 

0.97 

1.02 

1.08 

L.  Jia  et  al.  /  Journal  of  Power  Sources  253  (2014)  138-142 


141 


^  Sr(Sm)  Mn 

•  O  O  O  vacancy 

Fig.  3.  Top  view  of  the  oxygen  ion  diffusion  on  SSM  (100)  surface. 


atoms,  as  well  as  the  3d  states  of  the  relative  Mn  atom  before  and 
after  adsorption  are  shown  in  Fig.  4.  For  the  surface  without  O2 
adsorption  (see  Fig.  4(a)),  it  is  found  that  Mn  3d  orbitals  exhibits 
broad  bands  very  close  to  the  Fermi  level,  indicating  a  strong  itin¬ 
erancy  of  Mn  3d  electrons  [29],  thus  the  Mn4+  ions  can  easily 
transform  into  Mn3+.  There  is  an  overlap  between  Mn  3d  orbitals 
and  the  2 p  orbitals  of  the  adjacent  surface  O  atom  (03)  in  the  range 
of  -6  to  3  eV.  This  is  in  contrast  with  the  DOS  for  other  typical 
perovskites  with  large  band  gap  [25  .  And  the  overlap  reveals  that 
Mn  and  O  atoms  are  held  together  by  Mn-0  chemical  bond.  After 
adsorption,  as  shown  in  Fig.  4(b),  Mn  3d  electrons  transfer  to  lower 
energy  level  and  hybridized  with  the  2 p  electrons  of  the  adsorbed  O 
atoms.  The  hybridization  between  0  2 p  states  and  Mn  3d  states  is 
quite  strong  around  the  energy  level  of  -1.0  eV.  The  splitting  and 
broadening  of  peaks  indicate  a  strong  interaction  between  Mn  and 
the  adsorbed  O  atom,  which  facilitate  the  adsorption  of  O2  on  the 
atom  Mn  of  SSM50  surface.  The  overlap  between  the  2 p  states  of  01 
and  02  at  around  -1.0  eV,  -6.8  eV  and  -7.9  eV  reveals  that  01  and 
02  are  still  bonded  together.  It  can  be  seen  from  Fig.  2  that  the 
adsorbed  O  atom,  which  is  far  from  the  Mn  atom  (02)  has  an 
inclination  to  bond  together  with  03.  Thus,  the  2 p  states  of  02  and 


Energy  (eV) 


Fig.  4.  Calculated  PDOS  of  Mn  and  0  before  and  after  02  adsorption  on  SSM50  (100) 
surface.  The  Fermi  level  is  set  to  be  0. 


03  are  also  sketched  in  Fig.  4(b).  There  is  a  peak  formed  by  the  2 p 
states  of  03  at  about  -6.6  eV.  This  can  be  attributed  to  the  inter¬ 
action  between  02  and  03.  The  formation  of  02-03  bond  con¬ 
tributes  much  to  the  dissociation  of  the  adsorbed  O2  molecule. 

4.  Conclusions 

In  summary,  the  first-principles  density  functional  theory  cal¬ 
culations  were  performed  to  investigate  the  SSM  cathodic  reaction 
process:  bulk  diffusion  and  surface  adsorption.  It  can  be  concluded 
that  SSM50  is  a  potential  cathode  material  applying  for  IT-SOFC. 
The  obtained  results  are  as  follows: 

(1)  The  oxygen- vacancy  formation  energy  of  SSM  depends 
strongly  on  the  dopant  concentration.  The  formation  of  an 
oxygen  vacancy  in  SSM50  needs  an  energy  of  2.58  eV;  it  is 
easier  than  the  formation  in  LSM50. 

(2)  The  oxygen  migration  barrier  in  bulk  SSM50  is  0.75  eV.  The 
corresponding  diffusion  activation  energy  shows  that  SSM50 
is  the  most  suitable  cathode  material  among  these  systems 
considered. 

(3)  The  O2  molecular  can  be  adsorbed  on  the  SSM50  (100) 
surface  easily,  with  an  adsorption  energy  of  -3.94  eV.  The 
predicted  bond  length  of  O2  molecule  after  adsorption  in¬ 
dicates  that  the  absorbed  02  species  can  be  readily  disso¬ 
ciated  and  the  dissociation  energy  is  determined  to  be 
0.53  eV. 
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